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The one-dimensional discrete ordinates code ANISN has been adapted to 
transport low-energy (a. a few MeV) electrons. Calculated results obtained 
with ANISN are compared with experimental data for the transmitted electron 
energy and angular distributions for 1-MeV electrons normally incident on 
aluminum slabs of several thicknesses. The calculated and experimental 
results to date are in good agreement for a thin slab (0.2 of the electron 
range) but are not in good agreement for thicker slabs (0.6 of the electron 
range). Calculated results obtained with ANISN are also compared with cal- 
culated results obtained using Monte Carlo methods. 


INTRODUCTION 

The transport of low-energy (of the order of a 
few MeV) electrons through matter is important in 
the shielding of manned space vehicles that pass 
through the Van Allen electron belt. A code that 
treats this transport by means of Monte Carlo meth- 
ods is available (ref. l), but because of the poor 
statistical accuracy which can be obtained in some 
cases, a nonstatistical method of calculation is 
needed. To fill this need, the method of discrete 
ordinates has been adapted to the transport of low- 
energy electrons. The method of calculation is 
discussed in the next section, and the calculated 
results are presented and compared with experi- 
mental results and other calculated data in the 
last section. 

METHOD OF CALCULATION 


method has shown only partial success in transport- 
ing electrons. 

In the Monte Carlo treatment of electron trans- 
port, the individual electronic collisions are not 
considered, but rather the theories of multiple 
Coulomb scattering and continuous slowing down are 
used to group together large numbers of collisions 
(ref. 1). In the work reported here, the individ- 
ual electronic collisions are treated except that 
those electron-electron collisions which result in 
very small energy transfers (of the order of the 
average ionisation potential of the atom) are 
treated using the continuous slowing-down approxi- 
mation. The equation solved is the Boltzmann trans- 
port equation with a continuous slowing-down term 
added: 

$-V<|>(t,E,$) = p(3,e,s$) 


In principle, the discrete ordinates code 
ANISN (ref. 2) may be used to transport electrons 
by the simple expedient of introducing into the 
code the differential cross sections for electron- 
nucleus elastic collisions, electron-nucleus brems- 
strahlung-producing collisions, and electron- 
electron collisions. In practice, however, these 
cross sections are quite different from those which 
occur in neutron transport where the method of dis- 
crete ordinates has been used extensively , and the 

lThis work was funded by the National Aeronautics 
and Space Administration, Order H-38280A, under 
Union Carbide Corporations contract with the 
U. S. Atomic Energy Commission. 


E d 2 a (E’,E,fl'-S) 

+ n / ° dE- / dfl* 

E 


E 

/ c 

E+I ' 


+ J~° / dsr ln 


d 2 a. (E’,E,S’-n) 


dEdfi 


(S.E'.fr) 


n o T (E) 4>(^,E,n) 


+ (S(E) <(>($,E,3)] , 


( 1 ) 


816 



d2< V>e (E ' ’ E ’^' d2 ° el ( E ' . E »«' 4) 

dEdft dEdfi 

d 2 a tr (E',E,n’-2) 
+ dEdjj 


a e _^(E’) =! dE J dQ 


oT (E') = / dE / an 


d 2 o e+e ( E ' 4) 

dEdH 

d 2 o (E' .E.ft' 4) 

dEdfi 


o T (E') = o^(E-) + aJ n (E') , 


where 


R = a vector denoting the position 
of the particle; 

2 = a unit vector in the direction 
of the momentum vector ; 
dft = an element of solid angle ; 

E = the kinetic energy of an electron; 
= the electron flux per unit energy; 
E q = the highest kinetic energy con- 
sidered; 

n = the atomic number density; 

P(R,E,ft) = the number of electrons per unit 
energy per steradian per unit 
volume input at R from an ex- 
ternal source; 

the differential cross section 
for an electron (E f ) going in 
direction ft* to undergo elastic 
nuclear scattering (Coulomb scat- 
tering) to (E,ft); 
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the differential cross section 
for an electron (E T ,^*) to under- 
go bremsstrahlung (radiative) 
scattering with both nuclei and 
atomic electrons to (E,ft); 


dEdft 


the differential cross section 
for an electron (E * 1 ) to under- 
go inelastic scattering from 
atomic electrons to (E,i^), there- 
by producing a secondary electron 
with energy E'-E; 

I* = an arbitrary value taken to be 
the minimum energy loss allowed 


in an inelastic electron- 
electron collision; 

S(E) = the stopping power (energy loss 
per unit path length) due to in- 
elastic election-electron col- 
lisions resulting in an energy 
loss smaller than I* (these col- 
lisions are presumed to be 
straightahead) . 

The differential cross sections for electron- 
nucleus elastic collisions and electron-nucleus or 
electron-electron bremsstrahlung production are 
taken from standard sources (refs. 3 to 5). The 
differential cross section given by Moller (refs. 

3, 6) is used to describe electron-electron colli- 
sions which result in large energy transfers (> I*), 
The stopping power used to treat electron-electron 
collisions resulting in small energy transfers 
(< I 1 ) is obtained in the manner outlined by 
Rohrlich and Carlson (ref. 7) with the density ef- 
fect correction taken from Sternheimer (refs. 8, 9). 
The I' value used in the cases reported here is 
0.0163 MeV, i.e., 100 x I where I is the average 
ionization potential for aluminum, 163 eV (ref. 8). 

RESULTS AND DISCUSSION 

Goudsmit and Saunders on have obtained an ana- 
lytic expression for the angular distribution of 
transmitted electrons when monoenergetic electrons 
are normally incident on sufficiently thin slabs so 
that the energy degradation of the electrons may be 
neglected (ref. 10 ). Using the Goudsmit-Saunderson 
theory, Berger obtained the angular distribution of 
the transmitted electrons resulting from 1-MeV elec- 
trons normally incident on an aluminum slab of thick- 
ness 0.0287 g/cm 2 (ref. 11). The transmitted angu- 
lar current of electrons calculated by Berger is 
represented by the histogram shown in figure 1. Also 
shown in the figure are the results obtained for 
this same case with the discrete ordinates code ANISN. 
In the ANISN calculation, one energy group with a 
range from 1.0106 MeV to 0.989^ MeV and a midpoint 
of 1.0 MeV were used. No energy degradation was 
allowed and only elastic scattering was permitted; 
i.e., in this calculation equation 1 with a. = 0} = 

o/n\ 111 " br 

5(E; - 0 was solved. The two calculations shown in 
figure 1 are in excellent agreement, and thus the 
method of discrete ordinates can handle satisfac- 
torily small-angle multiple Coulomb scattering. 
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FIGURE 1.— Angular distribution of transmitted 
electron current for 1-MeV electrons normally 
incident on a 0 .0287-g/cm 2 -thick aluminum slab. 

Rester and Rainwater (ref. 12) and Rester and 
Dance (ref. 13) have measured the transmitted elec- 
tron current per unit energy from monoenergetic 
electrons normally incident on’ aluminum slabs of 
varying thicknesses. The solid histograms in fig- 
ures 2, 3, and 4 show the measured transmitted 
electron current per unit energy for aluminum slabs 
of thicknesses of 0.11 g/cm 2 , 0.22 g/cm 2 , and 
0.33 g/cm 2 , respectively. Also shown in the fig- 
ures are calculated results obtained with ANISN 
(solid curves) and calculated results obtained with 
the Monte Carlo transport code ETRA N of Berger and 
Seltzer (ref. l). In figure 2, the agreement be- 
tween the ANISN results and the experimental data 
is very good except at the very highest energies, 
and in this case, i.e., for a 0 .11- g/cm 2 -thick slab, 
the ANISN results are in slightly better agreement 
with the experimental data than the results given 
by ETRAN. In figure 3, neither ANISN nor ETRAN is 
in good agreement with the experimental data. In 
figure 4, which shows results for a relatively 
thick slab, 0.33 g/cm 2 (0.6 of the electron range), 
the Monte Carlo calculation is in very good agree- 
ment with the experimental data, while the ANISN 


calculation is in very poor agreement. The reason 
for the progressive failure of the continuous 
slowing— down ANISN calculation as. the thickness in- 
creases is not known. 
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FIGURE 2.— Transmitted electron current per unit 
energy per incident electron for 1-MeV electrons 
normally incident on a 0.11-g/cm 2 -thick aluminum 
slab. 
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FIGURE 3.— Transmitted electron current per unit 
energy per incident electron for 1-MeV electrons 
normally incident on a 0.22-g/cm 2 -thick aluminum 
slab. 
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FIGURE t .— Transmitted electron current per unit 
energy per incident electron for 1-MeV electrons 
normally incident on a 0 . 33-g/cm 2 -thick aluminum 
slab. 

The calculated and measured transmitted elec- 
tron current per unit energy and per unit solid 
angle is shown in figure 5 as a function of energy 
at several angles for a 0.11-g/cm 2 -thick aluminum 
slab. The calculated results agree fairly well 
with the experimental data at all angles, thus in- 
dicating an acceptable energy-angle correlation. 

Because of the lack of experimental data, it 
was not possible to present a comparison between 
calculation and experiment for the case of an elec- 
tron spectrum incident on a slab. Using the Monte 
Carlo code ETRM of Berger and Seltzer (ref. l), 
Scott (ref. Ik) calculated the transmitted electron 
current per unit energy for the case of a specific 
electron energy spectrum normally incident on alu- 
minum slabs, and this theoretical calculation was 
compared with results obtained with ANISN. The 
incident electron energy distribution used in the 
calculations is a representation of the spectrum 
resulting from thermal neutron capture in 23 5 U 

(ref. 15). This spectrum extends to electron 
energies of the order of 10 MeV and is shown ex- 
plicitly in reference Ik . 



E (MeV) 


-FIGURE 5 Energy distribution of transmitted elec- 
tron current at the specific angles indicated for 
1-MeV electrons normally incident on a 0.11-g/cm 2 - 
thick aluminum slab. 

The ETRAN results for an aluminum slab 0.5 g/ 
cm 2 thick are shown in figure 6 as a histogram and 
the ANISN results are shown as plotted points. A 
similar comparison is given in figure 7 for an 
aluminum slab 1.0 g/cm 2 thick. In both figures 6 
and 7 the ETRAN and ANISN results are in very good 
agreement. Thus, ANISN is apparently capable of 
treating accurately the transport of low-energy 
(< 10 MeV) electrons through relatively thick 
shields in some cases, e.g., the cases considered 
in figures 6 and 7> but not in all cases, e.g., 
figures 3 and 1*. 
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FIGURE 6 - Transmitted electron current P« ■ 

energy per incident electron for a specific energy 
spectrum (ref. 15) normally incident on a 0.5- 
g/cm z -thick aluminum slat. 
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FIGURE 7 .— Transmitted electron current per unit 
energy per incident electron for a specific 
energy spectrum (ref. 15) normally incident on a 
1 . 0 -g/cm z -thick aluminum slat. 
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